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1. Introduction 

A flat band is a dispersion relation that does not disperse. The kinetic energy of an 
excitation is suppressed, resulting in a Hamiltonian that admits a large degeneracy 
of localised, single particle eigenstates. The simplest examples of fiat bands are very 
familiar: Landau levels. However, in the absence of a magnetic field, they are more 
exotic beasts. Nonetheless, fiat bands can be manufactured in lattice models, typically 
by arranging a form of geometric frustration in which different hopping terms interfere. 

When flat bands are partially filled, interesting things happen. Arbitrarily small 
interactions or perturbations will lift the degeneracy, often driving the system to a 
highly correlated ground state. This mechanism provides the springboard for the rich 
spectrum of quantum Hall phases emerging from Landau levels. Moreover, in various 
lattice models, the flat band degeneracy has been used to drive instabilities towards 
itinerant ferromagnetism [1, 2] and superconductivity [3], as well as to fractional quan- 
tum Hall states [5, 6, 7, 8]. Flat bands have also made an appearance in studies of 
optical lattices [9], bilayer graphene [10] and topological phases [11]. 
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The purpose of this paper is to describe the emergence of a flat, dispersionless spec- 
trum for fermions in a holographic framework. Of course, in the context of AdS/CFT, 
highly correlated ground states are our bread and butter and we have little need to build 
them using flat bands. Nonetheless, we hope that the construction of a holographic 
flat band will allow a further, controlled study of this phenomenon. 

Concretely, we construct a non-relativistic boundary theory in d = 2 + 1 dimensions. 
However, in contrast to other gauge/gravity implementations of non-relativistic field 
theories, we do not change the bulk geometry. Instead we work in a bulk AdS4 spacetime 
and break Lorentz invariance only through the boundary conditions imposed on a Dirac 
spinor field. As we will see, it is possible to impose boundary conditions which preserve 
rotational invariance and scale invariance, but break both boosts and parity. The 
resulting 3d boundary theory can be thought of as a relativistic CFT, described by 
the 'W^ " degrees of freedom of the bulk geometry, coupled to a 3d Dirac spinor in a 
manner that breaks Lorentz invariance but, at least to leading order in 1/N , preserves 
z — 1 scale invariance. 

We shall see that the effect of the Lorentz breaking boundary conditions is rather 
brutal. The speed of propagation of the boundary fermion is driven to zero, resulting 
in an infinite fiat band, comprising of dispersionless, gapless excitations for all values 
of the spatial momenta. 

In Section 2 wc describe the novel boundary conditions, together with some features 
of the dual boundary field theory. Notably, the theory has a marginal, double trace 
operator, independent of the bulk fermion mass m. In some sense, this marginal op- 
erator can be thought of as the Goldstone mode for the broken Lorentz symmetry. In 
Section 3 we describe the single particle Green's function in vacuum and demonstrate 
the existence of the flat band. 

In Section 4, we place our theory at finite temperature and chemical potential by 
replacing the bulk with an AdS black hole. There has been much recent activity in the 
study of spinor two-point functions in these backgrounds due to the emergence of a 
Fermi surface revealing a rich spectrum of non- Fermi liquid behaviour [12, 13, 14, 15]. 
We find these Fermi surfaces numerically with our novel boundary conditions. We 
will see that, perhaps surprisingly, despite our singular spectrum, the formation and 
excitations of the Fermi surface are not greatly changed from the Lorentz invariant 
theory. 
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2. Non-Relativistic Fixed Points 



Gravity duals for non-relativistic field theories are typically implemented through the 
construction of a dual geometry that differs from AdS. The familiar examples are 
the Lifshitz geometry [18] and Schrodingcr geometry [19, 20], both of which describe 
quantum critical points with dynamical exponents z ^ 1. In this section we describe 
a non-relativistic system with z — 1 scaling symmetry in which Lorentz symmetry is 
broken by boundary conditions. 

To this end, we study a four-component Dirac fermion in asymptotically AdS4 back- 
grounds. We start in this section by considering vacuum AdS space of unit radius with 
Lorentzian metric 

ds^ — — ^ + r^rjui^dx'^dx" 

A detailed description of fermions in AdS spaces was presented in [21] whose conventions 
we (mostly) follow. The dynamics of a Dirac spinor is governed by the bulk action 



'bulk 



m 



(2.1) 



The covariant derivative is Dm = Om + \^ab,M^ where 0Jah,M is the bulk spin con- 



nection and r'*^ = |[r",r^]. In all these expressions, the capital index M denotes 
bulk spacetime while a, 6 denote bulk tangent space. The action (2.1) is not complete 
without boundary terms and we will spend some time shortly detailing the possibilities. 

To write the Dirac equation, it is convenient to choose the basis of bulk gamma 
matrices, 

r'^= M, A^ = o,i,2 , r = 
) 

where 7^ furnish a representation of the 3d Clifford algebra. We will later choose 
the basis 7^^ = (icr^, cr^, cr^). The Dirac equation can then be expressed in terms two- 
component spinors -0+ and each of which is an eigenvector of F^. 
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Working in Fourier modes for the boundary directions, ilj{x,r) = ■ilj{k]r)e^^'^ , the bulk 
Dirac equation is 



r ( rdr + - — m ] — —ik • 7 ■0. 



r ^r9r + ^ + "0- — ik ■ j ip+ 

Before discussing boundary conditions, it is useful to first examine the behaviour of the 
solutions near the AdS boundary r — > where, to leading order, 

i^-ik; r) = ^(fcjr-^/^-m ^ ^ _ _ _ as r ^ 00 (2.2) 

Here the 2-component spinors A, B, C and D are related by 

-'^(2^^1)3 , (2.3) 

We must decide which of these spinors is interpreted as the source and which is the 
response. As we now review, this is determined by the choice of boundary terms. 



2.1 Boundary Terms 

The bulk action (2.1) alone does not provide a full description of the dynamics. The 
problem lies with the terms which arise after integrating by parts. These are evaluated 
on the boundary of AdS 

^'S'buik = Q / d?'x\f—h (^6ijj+ijj- — 6ijj-ijj+ — ijj+Sijj- + ijj-Sijj+) + bulk term 
2 JdM 

Here h is the determinant of the induced boundary metric, h = gg^^ . 

Because the Dirac action is first order, it is not acceptable to insist that both S'^- 
and (5^+ vanish on the boundary. Restricting all components of the spinor in this way 
would be akin to fixing both position and momentum [22]. Instead we must fix just 
half of the spinor components which usually means either ip^ or ip^ although we shall 
discuss other options below. To achieve this, the bulk action must be augmented by 
a boundary term so that together they give rise to a well-defined variational principle. 
Depending on the value of the bulk mass m, there are a number of different ways to 
achieve this. 
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Standard Quantization 

For m > 0, the dominant term in (2.2) as we approacli tlie boundary is A{k). In tlie 
usual manner of liolograpliy, we would expect to fix the coefficient A{k), treating it as 
a source for the dual operator ^ in the boundary theory. This is achieved by adding 
the boundary action 

-^bdy = ^ / rf'xV^^V = ^ / d'xV^ + VS+V'-) (2.4) 

^ JdM ^ JdM 

The variation of the full action is then 

JdM 

which indeed vanishes if we impose Dirichlet boundary conditions for -0+. This choice of 
boundary condition is usually referred to as the standard quantization for fermions. The 
dual boundary field theory is a Lorentz invariant conformal field theory and contains 
a 3d Dirac fermion operator ^ , with dimension 




Alternative Quantization 

While the A{k) term in (2.2) is always dominant for m > 0, something rather special 
happens in the window < m < 1/2. Here the fall-off of A{k) becomes normalizable 
[21]. This means that while we are still at liberty to fix '0_|_ on the boundary if we want 
to, we are no longer obliged to do so since V'+ will happily fluctuate if we allow it. This 
provides us with an alternative choice of boundary condition in which we instead fix 
ip^ on the boundary. It is a simple matter to implement this: we need only add the 
boundary term with opposite sign, 

-^bdy = / d^X^/^jjjTjj (2.5) 

^ JdM 

which now results in a well defined variational principle if we impose Dirichlet conditions 
for '0_. This is known as alternative quantization. 

The ambiguity of boundary conditions for fields in AdS that lie within a particular 
mass window was first recognised in [23] and its implication for the AdS/CFT cor- 
respondence was explained in [24]: the two different boundary conditions correspond 
to two different dual CFTs on the boundary. In this alternative CFT, the fermion 
operator ^ has dimension 



Correlation functions in this alternative CFT are related to those in the standard 
quantization by a Legendre transform [24]. Note that as m — )■ 1/2, the dimension of 
the operator tends towards the unitarity bound A_ — > 1. Finally, note that taking the 
mass m < simply exchanges the role of the standard and alternative quantizations. 

Mixed Quantization 

If we are interested in Lorentz invariant boundary CFTs, then the two boundary condi- 
tions (2.4) and (2.5) exhaust the possibilities. Here we would like to relax this condition. 
Specifically, we will look for boundary conditions which preserve both the U{1) global 
symmetry tjj — )■ e^^ijj, rotational invariance and scale invariance, but not the full Lorentz 
invariance. 

There are two such choices, related by a discrete symmetry. The first can be imple- 
mented by the boundary condition 

•^bdy = 1 f d^xV^ r^r^V' I d^^^/^ Ulij+ + V'lV-) (2.7) 

^ JdM ^ JdM ^ ' 

With the basis of 3d gamma matrices — {ia^, a^, cr^), we write each two component 
spinor as 

The variation of the full action is now 

JdM ^ ' 

This vanishes if we impose Dirichlet conditions for V'+t and i)-^. 

The boundary field theory once again contains two fermionic operators, ^-f- and ^. 
The 3d rotation operator is i? = — 1[7^, 7^] = which means that ^'-|- has spin -|-| and 
has spin —\. However now the dimensions of these two operators differ, reflecting 
the broken Lorentz invariance. They are 

A[*t] = A+ , A[*J = A_ (2.9) 

Discrete Symmetries 

There is a second mixed quantization where the boundary condition (2.7) differs by an 
overall minus sign, requiring that we impose Dirichlet conditions for ■0+4, and i)-^. The 
resulting boundary operators have dimension A[^-^] = A_ and A[^j,] = A+. 
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These two conformal field theories are related by a parity transformation. Recall 
that in d = 2 + 1 dimensions, parity acts by reflecting just a single spatial coordinate, 
— > —x^ with the corresponding action on spinors P : '^^ ^ Thus, while the 
Lorentz invariant theories preserve parity, it is broken by our non-relativistic bound- 
ary conditions^. The non-relativistic boundary theory also breaks charge conjugation, 
which acts as C : ^ ^^1. 

2.2 A Poor Man's RG 

Dual CFTs described by different boundary conditions are not unrelated. There exists 
a renormalization group (RG) flow from the alternative quantization to the standard 
quantization, initiated by turning on a relevant double trace operator [25, 26]. 

Some properties of this RG flow for fermions were described in [27]. One starts in 
the alternative quantization and perturbs the boundary theory by the Lorentz invariant 
double trace operator '1'^. This is a relevant operator, with dimension A['5'^] = 3 — 2m. 
The end point of the flow is the standard quantization. For this reason, the alternative 
quantization is sometimes referred to as the UV CFT, while the standard quantization 
is the IR CFT. 

What about our non-relativistic dual CFT? This too can be reached by an RG flow 
from the alternative quantization. If we care only about rotational symmetry and 
the U{1) global symmetry, the UV CFT contains two relevant operators: and 
*7°*. To flow to the non-relativistic theory with operators (2.9), we want to change 
the boundary conditions for the bulk field leaving ip^^ untouched. This can be 
achieved by perturbing with the double trace operator 

J (fx + i7°)^ = j d^x 2^1^^ (2.10) 

It is useful to remind ourselves what these two operators would mean for a free fermion. 

is a mass term; it breaks parity, but respects charge conjugation. Meanwhile, 
ijf^o^ is a chemical potential for a free fermion (at least in a suitable thermodynamic 
ensemble). It breaks charge conjugation but respects parity. Had we added the operator 
(2.10) for a free fermion, it would open up a gap in the spectrum, with the chemical 
potential tuned to sit at the bottom of the band, resulting in a single component, non- 
relativistic, gapless fermion. Of course, our fermions are far from free. We will see 
shortly what becomes of them. 

^We note that we could equally well implement these mixed boundary conditions for spinors in 
Lifshitz geometry. Since there is no Lorentz symmetry of the boundary theory in this case, the only 
price paid is the breaking of parity. 
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Figure 1: Flows between the UV, IR and non-relativistic fixed points. The red circle denotes 
a of fixed points, breaking the U{1) global symmetry. Parity reflects the diagram about 
the vertical. 

Perturbing the UV CFT by ^'j^E'i- will result in a flow to the other non-relativistic 
flxed point. We leave a full analysis of the flow equations to future work [28]. For 
now, it is comforting to know that we can reach the non-relativistic flxed point from a 
sensible starting point, albeit at the cost of flne tuning an RG flow. 

Each of the non-relativistic flxed points also contains a relevant operator, namely 
"^l^^l (or, in the parity related theory, In each case, turning on this operator 

takes us back to the IR CFT associated to the standard quantization. The picture of 
two dimensional RG flows that emerges is shown in Figure 1. 

A Circle of Fixed Points 

There is one further double trace operator of interest at the non-relativistic flxed point, 
namely the charged, rotationally invariant operator \E'i-\I'4.. Usually the dimensions of 
operators are determined by the mass of a bulk fleld. Not so for this operator: its 
dimension is + A_ = 3. In other words, it is marginal for all values of the bulk 
mass in the window |m| < 1/2. 

A similar mechanism for generating a marginal operator in the case of complex scalar 
flelds was described in [25]: the real part of the scalar was given standard boundary 
conditions; the imaginary part alternative. In both cases, the marginal operator can 
be loosely thought of as a Goldstone mode for the symmetries broken by the boundary 
conditions. For us, the broken symmetry is Lorentz symmetry. 
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Turning on the marginal operator ^^-^j, sweeps out a manifold of fixed points; this 
is the red circle in Figure 1. (For readers without colour, it is the circle in Figure 
1). Since this operator is complex, there is a two-dimensional set of fixed points. The 
magnitude of the coefficient takes us around the circle; the phase (not shown in the 
figure) takes us from the back of the circle to the front. The result is that the manifold 
of fixed points is topologically S^. Since we are turning on a charged operator, each of 
these CFTs away from the two special points breaks the f/(l) global symmetry of the 
boundary theory. In fact, we can define these fixed points directly by introducing the 
linear combinations, 

= COS e V'it + sin , X±i = cos 6 V'±i - sin 6 
It is then simple to check that the addition of the boundary term, 

^bdy I d'x^ (xlx+ + Xlx-) (2.11) 

^ JdM ^ ^ 

requires us to impose Dirichlet conditions for x+-|- and X-\r- Because x+-|- is a linear 
combination of -i/^+i- and both of which have the same radial dependence, these 
boundary conditions preserve scale invariance. However, as promised, they break the 
global symmetry apart from at the poles, ^ = and Q — 7r/2, where these boundary 
conditions reduce to our previous mixed quantization (2.9) and its parity twin. 

Finally, the RG properties of a massless bulk fermion, m = 0, are worthy of comment. 
In this case, there is no RG fiow from alternative to standard quantization since both 
and i|f7°^ are marginal. In particular, the marginal operator gives rise to a 
circle of fixed points that interpolate between the alternative and standard quantization, 
as previously described in [29]. This is again a "Goldstone operator", associated to 
the breaking of bulk chiral symmetry by the boundary condition. Together with the 
Lorentz violating 'J'7°\E', the two operators sweep out a torus of charge invariant 
fixed points. Further including the charge violating operator ^1-^4,, results in a four- 
dimensional manifold of fixed points. 

3. Green's Functions 

We now turn to the excitation spectrum of the non-relativistic boundary theory. We will 
first provide a quick method to determine the Green's function. However, as previously 
advertised, the resulting spectrum will have the rather peculiar property of a fiat band 
and for this reason we subsequently provide a second derivation. 
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It will prove useful to first recap the fermion propagator for relativistic conformal 
theories. The two point function, which is fixed by symmetries, was first computed 
holographically in [30, 31], with various properties of the retarded propagator clarified 
in [21]. In this case, the Dirichlet boundary conditions for some components of the 
spinor at the boundary r — )■ oo are accompanied by ingoing boundary conditions at the 
AdS horizon r — 0. The Dirac equation is then easily solved in terms of appropriate 
Bessel functions. 

As we reviewed in Section 2, for the standard quantization ■0+ is fixed on the bound- 
ary, and the two-spinor A in (2.2) is interpreted as the source, while D is the response. 

The relation between the two is written D[k) = S{k)A{k). The retarded propagator is 
then defined to be G/j(A+) = —iS'j'^, where the factor of 7'^ in this expression reflects 
the fact that the Green's function computes (^^^) rather than (^f^). For timelike 
momenta, a; > |A;|, one finds the propagator 

Gr{A+) = aA+ {co' - \k\Y^-'{k ■ 7)7° (3.1) 

where the overall constant is given by 

_ e-(^-^)" r(2 - A) 
"A - 22A-3 r(A - 1) 

In the alternative quantization, we instead fix with D now interpreted as the source 
and A the response. To compute the propagator, it is not necessary re-solve the Dirac 
equations in the bulk. Wc need only realise that the roles of the source and response 
have been interchanged so that the propagator is given by ^^^(A-) = iS~^'y^. Indeed, 
as a check one can easily verify that the two-point function obeys 

3.1 Non- Relativistic Boundary Conditions 

With mixed, non-relativistic boundary conditions, the sources are A^ and D^. The 
propagator is again easily determined from (3.1) without going through the rigmarole of 
solving the bulk Dirac equations. We need only massage the equation D{k) = S{k)A{k) 
into the form. 




The retarded propagator is Qr — —S. The lack of Lorentz invariance in the theory 
means that there is no 7° factor in this relationship and it is simple to check that 
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this definition is consistent with the requirement of unitarity since it ensures that the 
imaginary part of the diagonal components of the propagator are positive. For timehke 

— * 

momenta, co > \k\, the retarded propagator is 



«A+(u;2-|^|2)^+-i iki-k2 



(3.2) 



For oj < —\k\, the coefficients q;a± in the above expression should be replaced by 
a\_^. It is perhaps unusual to see a two-point function mixing operators of different 
dimensions but this is not forbidden in a non-relativistic system. (The usual proof that 
this cannot happen makes crucial use of special conformal transformations). However, 
by far the most unusual aspect of this propagator is the 1/uj pole out-front. This is 
the first indication of a fiat band. This pole survives at spacelike momenta, > a;^, 
where the retarded Green's function is given by 



uj + ie\ -iki-k2 \aA.\i\k\'' - u'')^--^ 



(3.3) 



We have added the +ie prescription to ensure that the propagator is analytic in the 
upper-half complex ou plane, as required by causality. We will see later that this is 
consistent when we study our system at finite temperature. 

To better understand the structure of the propagator, let us look at the eigenvalues, 
which we will denote as A+ and A_. The propagator obeys det Qr — —1, which tells us 

— * 

that the eigenvalues come in reciprocal pairs. For small u <^ the two eigenvalues 
scale as A+ ~ —l/ou and A_ ~ ou. Specifically, 

X+^-^ \\aAAP"' + and A_ ^ ( |q;a+ l^^"^ + ^ 




|Q!A+|A;2"*y \k\ \ |q;a+|A;2'"^ 

For the special case of m = the eigenvalues simplify to X± = u^^Vk"^ — u"^ ± |A^|). 
The corresponding eigenvectors are {\k\,^{iki — k2))- 

One of the eigenvalues exhibits the peculiar property of a pole at a; = for all values 
of the spatial momentum k. This implies the existence of localised, gapless excitations 
with vanishing speed of propagation. As discussed in the introduction, dispersionless 
modes of this kind are referred to as a flat band although, in the absence of a magnetic 
field, it is more usual to find such a spectrum in lattice models where the band covers 
only a finite interval of momenta. The holographic theory exhibits an infinite fiat band. 
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The appearance of a flat band is very surprising. To illustrate how extreme this 
behaviour is, one only has to note that the very language of physics breaks down: 
dispersion relations do not disperse; propagators do not propagate. It is a form of local 
criticality, in which neighbouring spatial points no longer speak, but different from the 
2; — >■ oo, AdS2 criticahty that arises at finite density [15]. Needless to say, it would be 
interesting to understand the emergence of this flat spectrum from the perspective of the 
boundary theory. One clue may come from the analysis of a free fcrmion. As discussed 
in the previous section, the deformation (2.10) gaps the spectrum of a free fcrmion 
and leaves the chemical potential kissing the bottom of the band, resulting a gapless 
excitation with quadratic dispersion. For our strongly interacting fermion, such a 
quadratic dispersion is not possible since the coupling to the background AdS spacetime 
forces the correlation function to enjoy z — 1 scaling with an associated uj — vk 
dispersion. Thus we seem to be in a situation where the deformation (2.10) removes 
the linear dispersion term but higher order terms are forbidden by the background, 
driving the speed of propagation v — > 0. 

One might worry about the existence of flat bands in a theory which allows for 
pair creation. If we could populate the band with particle anti-particle pairs from 
the vacuum, then our theory would be rabidly unstable. Thankfully, that is not the 
case. The excitations of the band have fixed charge under the global symmetry of the 
boundary^ and charge conservation requires that any excitation of the flat band with 
eigenvalue A+ that is created from the vacuum must be accompanied by its partner 
with eigenvalue A_. Fortunately, this means that the flat band cannot be populated 
from the vacuum. We must place particles in by hand. 

As mentioned in the introduction, a flat band is most interesting because it is unstable 
to interactions and perturbations, often driving the physics to an interesting strongly 
coupled corner. In the next section, we will populate and heat our band and find that it 
is surprisingly stable, at least when the fermions are treated in the probe approximation. 
The discussion above suggests that the flat band is most likely a large N artifact and 
the most interesting question is how, if at all, it is lifted by 1/A^ corrections. One such 
source of corrections arises from quantum fermion loops in the bulk. In particular, 
since the boundary conditions break Lorentz invariance, bulk vacuum bubbles, which 
usually only renormalize the cosmological constant, are now at liberty to change the 
bulk metric away from AdS^. 

^Strictly, the boundary theory has a global symmetry only if we gauge the V' — >■ e'"V symmetry in 
the bulk. We implicitly assume this is the case. We will be more explicit in the following section. 
^We thank D.T. Son for discussions on this point. 
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Finally, we would like to mention that something similar to a fiat band was seen 
already in holographic electron star calculations, where the Fermi surface is smeared 
into a Fermi ball, with gapless excitations at all points in the interior [32] . This ball is 
resolved hy 1/N corrections into a family of discrete nested Fermi surfaces [33, 34]. 

3.2 Bulk Zero Modes 

The derivation of the Green's function presented above used the relationship with the 
relativistic theory. One can also extract the Green's functions in a more pedestrian 
manner directly from the bulk solutions. This method has the advantage of exhibiting 
the normalizable bulk modes corresponding to the a; = pole. 

The solution to the bulk Dirac equation is well known [31]. Restricting to spacelike 
momenta, /c^ = — u;^ > 0, the normalizable solution is given by the Bessel functions 

where ^± = ^±771. Expanding to leading order in r, the two-component spinors x+ 
and X- ^'I'c related to the spinors A and D in (2.2). However, we must still satisfy 
(2.3) which provides a relationship between source and response. In this manner, it is 
a simple matter to reproduce the Green's function (3.3). 

Poles in the Green's function correspond to normalizable bulk modes in the absence 
of a source. Thus, the fiat band Green's function (3.3) tells us that the mixed boundary 
conditions should allow for u = modes for arbitrary k. Turning off the source sets 
— = 0. The conditions (2.3) then requires that the ip^ and ip^ solutions are 
related by 

\ / uj ki-ik2\ I 
J \k1-ik2 -oj / \ 

When u = 0, this equation can indeed be solved for arbitrary k. One can also check 
that the bulk Hamiltonian, "H ~ iplY^r + ^i7* ~ fn)ijj = evaluated on these solutions, 
courtesy of the spinor contractions. 

4. Black Hole Backgrounds 

In this section we discuss the behaviour of our non-relativistic fermion in the presence 
of an AdS black hole, corresponding to finite temperature and finite chemical potential. 
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One may expect that as soon as the bulk theory is placed at finite temperature, all 
modes in the flat band will become occupied. In fact, we will see that, at least to 
leading order in 1/N, this is not the case. 

The AdS Reissner-Nordstrom black hole metric is 

where the emblackening factor is given by 

1 + Q\ 



We have chosen to rescale the coordinates to place the black hole horizon at r = 1 which 
means the temperature and chemical potential described below are both dimensionless. 
The background also contains an electric field, with gauge potential given by 

where the gauge coupling l/g"^ is the coefficient in front of the bulk Maxwell action. In 
the context of AdS/CFT, this gauge coupling arises in the current two-point function, 
{J J) ~ i/g^i where it is loosely interpreted as the amount of charged stuff hving 
among the N'^ strongly coupled degrees of freedom of the boundary field theory. In 
what follows, we will refer to this N"^ degrees of freedom as the soup. 

As usual, this black hole background places the strongly coupled soup at finite tem- 
perature T and chemical potential /i, given by 

We probe this soup with our fiat band fermion. The bulk action (2.1) remains the 
same, but the covariant derivative now conceals a coupling to the background gauge 
field Dm — Om + \(^ab,M^"'^ ~ iQ^M- We will set the charge q — I and study how the 
physics of the probe changes as wc vary the gauge coupling g. Note that, by a rescaling 
of the gauge field, this is entirely equivalent to varying q. However, it is somewhat 
simpler to interpret variations of the gauge coupling: as g increases, the amount of 
charged matter sitting in the soup decreases and the probe fermion gains prominence. 
In what follows, we work exclusively with a fermion of mass m = 0. For masses m ^ 0, 

— * 

it is less easy to separate the eigenvalues A+ and A_ for all values of k and u but 
invariant objects like TrQn do not exhibit qualitatively different behaviour from the 
m — story we tell below. 



14 



Figure 2: The imaginary (left) and real (right) parts of the flat band eigenvalue, A+, at 
finite temperature 

Finite Temperature 

We start by setting g = 0, corresponding to placing the boundary field theory at finite 
temperature, but vanishing chemical potential. The spectral function for the singular, 
flat band eigenvalue A+ is shown in the left hand plot of Figure 2. Strikingly, although 
there is a depletion in the low momentum modes, the spectral function for the high 
momentum modes is essentially unchanged from the vacuum. In particular, there is 
an w = peak which becomes strongly delta functionesque as k increases as evidenced 
by the pole in the real part of the Green's function which is plotted on the right of 
Figure 2. It is unusual to see such dramatic response at energies below the temperature 
and suggests that there are zero energy fermionic modes which have not been excited 
even though the system is at finite temperature. These modes persist because they sit 
outside the hghtcone and cannot therefore decay. 

To understand physically what is happening here, we need to realise that the black 
hole background does not directly heat up the fiat band. Rather, it heats up the 
CFT soup to a temperature which, in our convention, is T = 3/Att. This hot soup 
is subsequently probed by the bulk fermion. Usually, one expects that the probe and 
soup equilibrate, so that all probe states up to energy T are excited. Yet this does not 
happen for the fiat band, where u = states remain unoccupied. It would seem that 
the question of equilibrium is more subtle. Specifically, the soup is relativistic and has 
momentum modes excited up to \k\ ~ T. Since momentum is conserved, these soft 
modes can only excite the high momentum modes of the fiat band under the unlikely 
condition that many collide at once. Thus, the spectral function of the fiat band is 
determined not just by the temperature, but by the momentum distribution of the soup 
and its coupling to the fermion. 
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Figure 3: The imaginary parts of the eigenvalues A+ (left) and A_ (right) for g = 1 



Im[,\^i 




Figure 5: The imaginary parts of the eigenvalue A_ for g = 4: 
Finite Density 

We now turn on the chemical potential, focussing on T = 0. In Figure 3, we plot 
both eigenvalues of the spectral function, ImA± for g = 1. In the left-hand figure, 
ImA_|_ clearly shows the fiat band. Because frequency is now measured relative to the 
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Figure 6: The density plot of ImA-|_ (left) and ImA_ (right) as g varies. These plots were 
made at cj = — 10~^ and are robust to lowering the value of |t<;|. 

chemical potential, the band is shifted to w = — /i = —y/Sg. Once again, we see that 
strong peaks remain at high momenta, reflecting the fact that these states are not 
occupied. (The apparent disappearance of the peak at high k is numerical artifact. It 
arises because the peak tends sharper as can be checked by plotting a suitable cross- 
section at fixed k). This means that, despite equality in their energies, the flat band is 
occupied from the low momentum modes up and, for this reason, does not contribute to 
the ground state degeneracy in the background of the Reissner-Nordstrom black hole. 
For low momenta, the flat band becomes mildly dispersive, albeit with a large width. 
For g = 1, the spectral function for A_ is essentially featureless. 

In Figures 4 and 5 we crank up the gauge coupling to g = 4 and, to illustrate the 
features, present a birds eye view of the spectral function together with the landscape 
plots. Once again, the flat band is clearly visible in Figure 4. However, now both 
eigenvalues A+ and A_ reveal a Fermi surface. The presence of the flat band does 
not seem to greatly alter the Fermi surface. Indeed, it can be checked that the quasi- 
particle excitation spectrum (or lack thereof) is governed by the AdS2 near horizon 
regime of the black hole, as explained in [15]. At most, the flat band appears to mildly 
suppress the formation of the Fermi surface. For example, note that in both standard 
and alternative quantizations (which are effectively the same for m = 0), the Fermi 
surfaces is already present at g = 1. Moreover, as we increase g, the Fermi surface 
appears in the A_ eigenvalue before it appears in A+. 

This latter point is highlighted in Figure 6 which depicts the a; = density plot as 
g is varied. Above the black, diagonal line lies the oscillatory region, defined in [15], 
where the IR exponent is complex and no Fermi surface can form. The Fermi surfaces 
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Figure 7: The quasi-normal mode spectrum for 5 = 4 and A: = 0.1, A; = 0.3, k = 5 and k = 7 

appear when the white hnes extend below the black, where they are delta functions 
revealing themselves as faint, numerical pinpricks of paint. Note also the splash of 
white along the g = axis in the left-hand plot, corresponding to the flat band. 

Finally, in Figure 7 we plot the quasi- normal mode spectrum for T = 0, (7 = 4 and 
increasing values of the momentum. The plots show k = 0.1, k = 0.3, k = 5 and k = 7 
respectively. They were computed using an extension of the scalar method presented 
in [35]. (See also [36]). Here one Taylor expands the Dirac equation in the black 
hole background, transforming the differential equation into a matrix equation. This 
expansion is truncated at finite order, resulting in a large but finite matrix and, for fixed 
k, one scans the complex u plane for points where the determinant vanishes. These are 
the quasi-normal modes, plotted for both A+ and A_ in the diagrams above. (These 
diagrams were made after truncating to a 200 x 200 matrix). As the quasi-normal 
modes approach the real axis, they coincide with peaks in the spectral function. 

The first feature in these plots is the ridge of quasi-normal modes running almost 
vertically downwards. These are the remnant of a branch cut. (Indeed, as we go to 
higher order, we more points appear in this ridge). For small \k\, a number of quasi- 
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normal modes fan out from this ridge in the Recu < direction. These are clearly 

— * 

visible in the first and second plots above. As |A;| increases yet further, a number of 
these quasi-normal modes start to approach the real a; axis. For the g — 4: case shown, 
the first two quasi-normal modes cross over to Rea; > where they hit the real axis 
before bouncing off. These correspond to the Fermi surfaces shown in 4 and 5. The 
first quasi-normal mode to hit is the A+ eigenvalue; the second is A_. Indeed, it can 
be verified that A+ exhibits a Fermi surface at lower kp which indeed exhibits a Fermi 
surface at lower \k\ than A_. The third quasi-normal mode does something different: 
it too heads towards the real u axis, but it doesn't bounce. It can be seen in the last 
two plots above, nudging closer to uj/Q = —g = —4 where it remains for ever higher k. 
This is the fiat band. 
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